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摘 要 : 疏 勒 河 流域 地 处 河西 走廊 最 西端 ,是 典型 的 干旱 内 陆 河 流域 ,生态 环境 极其 脆弱 ,也 是 我 国 
重要 的 生态 屏障 之 一 。 借 助 Fragstats 软件 从 斑 块 和 景观 2 个 尺度 ,结合 2000 一 2018 年 土地 利用 数 
据 ,评价 分 析 聋 勒 河 流域 景观 生态 风险 时 空 变化 特征 ,利用 地 理 探测 器 定量 分 析 了 该 流域 景观 生 
态 风险 驱动 因素 。 结 果 表 明 :(1) 疏 勒 河流 域 主要 景观 类 型 为 未 利用 地 和 和 草地 ,其 斑 块 个 数 (NP)、 
景观 形状 指数 (LSI) 、 最 大 斑 块 指数 (LPI) 和 聚集 度 指数 (AI) 值 均 较 高 ;2000 一 2018 年 蔓延 度 指数 
(CONTAG ) 值 下 降 、 香 农 多 样 性 指数 (SHDI) 和 香农 均匀 度 指数 (SHEI) 值 缓慢 上 升 ,流域 景观 破碎 
度 严 重 。(2) 朴 勒 河流 域 景 观 生态 风险 空间 分 布 呈 现 北 高 南 低 ,2000 一 2018 年 呈现 流域 景观 生态 风 


险 逐 渐 下 降 趋势 ,其 中 较 高 风险 和 高 风险 地 区 了 


殿下 降 明 显 。(3) 人 为 干扰 度 是 影响 景观 生态 风险 


空间 分 布 的 主要 因子 , 归 一 化 植被 指数 (NDVI) 次 之 ,人 口 密 度 影 响 最 小 ;景观 生态 风险 的 影响 均 为 


双 因 子 增强 或 非 线性 增强 ,无 显著 差异 大 多 表现 在 


自然 和 自然 因子 间 ,自然 和 人 为 因子 之 间 存 在 


显著 性 差异 。 因 此 ,对 流域 进行 景观 生态 风险 指数 评价 和 驱动 因素 分 析 至 关 重 要 。 
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景观 生态 风险 是 景观 格局 和 生态 过 程 相 互 作 
用 产生 不 利 的 影响 ,主要 表现 在 人 类 活动 或 自然 灾 
害 对 景观 的 组 成 .结构 功能 等 影响 。 流 域 是 连 
接地 表 水 .土地 植被 覆盖 和 自然 -社会 的 重要 综合 
生态 地 理 区 域 ,也 是 最 为 重要 的 地 理 景观 , 是 受 
到 人 类 活动 强烈 作用 的 区 域 ,流域 保护 是 实现 区 域 
社会 经 济 可 持续 发 展 的 必然 选择 ”。 自 然 因素 和 人 
类 活动 已 经 对 景观 格局 产生 了 深远 的 影响 ,流域 
尺度 上 的 景观 生态 风险 分 析 评 价 是 对 流域 进行 生 
态 环 境 保护 的 重要 研究 内 容 , 也 成 为 生态 风险 评 
价 的 热点 之 一 ”。 国 内 外 学 者 分 别 在 国家 *”、 区 
JU 乡镇 "流域 "等 不 同 尺度 上 对 景观 生态 
风险 评价 的 理论 与 方法 进行 了 广泛 探讨 ,已 初步 形 
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成 了 研究 框架 "。Yang 等 一 在 改进 TOPSIS 模型 基 
础 上 ,选取 了 24 个 指标 从 农业 .工业 和 城市 化 3 个 
方面 分 析 评 价 渭河 流域 生态 风险 ,认为 通过 改善 水 
质 和 保持 生物 完整 性 可 以 调控 该 流域 生态 风险 。 
李 青 半 等 5 采用 最 小 积累 阻力 模型 和 空间 主 成 分 
分 析 对 宁 江 进行 景观 生态 风险 评价 ,研究 发 现 人 类 
社会 活动 和 景观 格局 对 景观 生态 风险 影响 较 大 , 模 
拟 发 现 通过 优化 景观 格局 使 得 景观 连通 度 显 著 增 
加 ,可 以 改善 景观 生态 风险 程度 。 郑 杰 等 史 利 用 遥 
感 影像 研究 2009 一 2017 年 草 海 流域 景观 生态 风险 
演变 特征 ,结果 表明 在 草 海流 域 高 生态 风险 、 较 高 
生态 风险 和 中 生态 风险 面积 均 呈 不 同 程度 增加 态 
势 , 旦 景观 生态 风险 存在 显著 的 空间 自 相 关 性 。 整 
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孙 丽 区 等 : 基于 地 理 探测 器 模型 的 下 勒 河流 域 景观 生态 风险 评价 及 驱动 因素 分 析 


体 来 看 , 当前 所 用 研究 的 景观 因子 大 多 集中 在 对 景 
观 格局 的 量化 和 分 析 ,利用 研究 结果 直接 或 间接 的 
反映 景观 生态 风险 。 对 景观 生态 风险 驱动 力 研 究 
多 从 定性 方面 入手 ,缺乏 定量 研究 手段 。 地 理 探 
i d 77 是 一 种 基于 统计 学 原理 ,通过 分 析 层 内 方 
差 与 层 间 方 差 的 差异 定量 表达 人 研究 对 象 空间 分 异 
性 的 研究 方法 ,该 方法 不 仅 可 以 探测 数据 型 数据 和 
类 型 数据 ,也 可 以 定量 探测 驱动 因子 之 间 交 互 作 用 
以 及 判断 每 2 个 因子 对 景观 生态 风险 空间 分 布 的 影 
响 是 否 具 有 显著 差异 ,与 其 他 空间 异 质 性 探测 工 
有 具 相 比 地 理 探 测 器 具有 更 高 的 解释 效率 。 目 前 地 
理 探测 器 广泛 应 用 于 定量 研究 气候 变化 ”” CO, 
减 排 “、 植 被 变化 “居民 点 分 布 " 力 至 新 冠 
(COVID-19) 疫 情 散布 ”的 驱动 因素 ,然而 在 分 析 
景观 生态 风险 驱动 因子 中 尚 不 多 见 , s SCELTE ELA] 
河流 域 景 观 生 态 风 险 评 价 基础 上 ,利用 地 理 探测 器 
对 其 景观 生态 风险 驱动 因素 开展 定量 分 析 。 

草 勒 河流 域 是 典型 干旱 内 陆 河流 域 ,生态 环境 
脆弱 ,近年 来 受到 气候 变化 和 人 类 活动 影响 , 琉 勤 
河流 域 景观 趋 于 破碎 化 ” ,导致 景观 生态 风险 进 一 
步 加 剧 ,深入 研究 该 流域 景观 生态 风险 时 空 演变 特 
征 , 识 别 不 同 驱 动因 子 对 景观 生态 风险 的 贡献 量 ， 
以 期 为 控制 琉 勒 河流 域 景观 生态 风险 升级 ,以 及 流 
域 优 化 管理 和 协调 可 持续 发 展 提供 科学 依据 。 


1 研究 区 概况 


玻 勒 河流 域 位 于 我 国 西北 干旱 区 ,是 我 国 唯一 
一 条 由 东 向 西 的 内 陆 河 ,处 于 河西 走廊 最 西端 , 介 
F 92?22' .98935'E,38*05' -42*45' N 之 间 ,流域 总 面 
fR25 1.12x10? km2:( 图 1)。 路 勒 河 发 源 于 祁连山 脉 
的 岗 格 尔 肖 合力 岭 , 流 经 肃 北 蒙 古 族 自治 县 .敦煌 
市 、 玉 门市 以 及 瓜 州 县 。 流 域 地 势 北 低 南 高 ,海拔 
介 于 802~5504 m。 流 域内 降水 量 少 , 蒸 发 量 大 , 东 
西 2 侧 临近 沙漠 ,水 源 补 给 主要 为 冰川 融 水 和 降 
水 。 上 游 地 区 海拔 高 ,地势 陡 峭 , 存 在 永久 性 冰川 
雪 地 ,中 下 游 地 势 平 坦 ,绿洲 与 荒漠 并 存 , 其 中 未 利 
用 地 面积 占 70% 以 上 ,土地 沙漠 化 问题 严重 。 根 据 
《2019 年 甘肃 省 水 资源 公报 》,2019 年 芍 勒 河流 域 有 
人 口 5.36x10; 人 ,GDP 为 3.67x10" 元 ,耕地 面积 》 
14.79x10* hm。 
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Fig. 1 Spatiallocation of the Shule River Basin 


2 数据 与 方法 


2.1 数据 来 源 及 处 理 

2000.,2010 2018 年 土地 利用 数据 以 及 县 级 行 
政 区 划 矢 量 数据 来 源 于 中 国 科 学 院 地 理科 学 与 资 
源 研究 所 资源 环境 数据 云 平台 (http://www.resdc.cn/ 
Default.aspx) ,空间 分 辨 率 为 1 km ,投影 坐 标 系 均 为 
Krasovsky_1940_Albers。 其 中 ,耕地 城市 及 农村 居 
名 点 的 数据 采集 ,正确 率 不 低 于 95% ,草地 、 林 地 、 
水 域 等 正确 率 不 低 于 90% ,未 利用 土地 正确 率 不 低 
于 85%。 参 考 国家 土地 利用 现状 分 类 标准 ,适当 合 
并 后 确定 耕地 、 林 地 、 草 地 水域 .居民 地 和 未 利用 
土地 6 类 为 本 次 研究 的 景观 类 型 。 选 取 的 驱动 因子 
为 3 个 气候 因素 (年 降水 量 、 年 最 高 气温 \ 年 最 低 气 
温 )、3 个 地 表 因 素 [ 高 程 (DEM) .坡度 . 归 一 化 植被 
指数 (NDVI) ] 和 2 个 人 为 因子 (人 口 密 度 和 人 为 干 
扰 度 )。 其 中 人 口 密度 来 源 于 Worldpop 网 站 (http:// 
www.worldpop.org/ ) ;年 降水 量 、 气 温 数 据 来 源 于 国家 
地 球 系统 科学 数据 中 心 ;NDVI 和 DEM 来源 于 中 国 
科学 院 地 理科 学 与 资源 研究 所 资源 环境 数据 云 平 
f ,坡度 数据 则 由 DEM 数 据 在 GIS 中 提取 。 

年 降水 量 ,高程 NDVIL 气温 、 人 为 干扰 度 在 
ArcGIS 中 用 自然 断 点 法 分 为 10 类 ,人 口 密度 用 等 数 
量 间距 法 分 为 10 类 ,坡度 按 术 土 保持 综合 治理 规划 
通则 (GB/T 15772 一 2008) 分 为 6 类 (0~3°、3°~8° 、8°~ 
15° .15°~25° .25°~35° 和 >35°)。 
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22 研究 方法 

2.2.1 景观 格局 指数 ”景观 格局 指数 能 够 从 无 序 的 
景观 斑 块 中 发 现 潜在 的 规律 及 意义 , 即 用 定量 化 的 
指标 来 反映 景观 斑 块 在 景观 组 成 结构 和 空间 上 的 
RIE, RSZ ENEA RRA A 
EREE 8 4) B AR DU, MARME RK F 
上 针对 个 数 、 边 长 和 分 布 聚 集 性 和 连接 性 多样 性 
等 方面 ,利用 Fragstats 选取 斑 块 个 数 (NP) .最 大 斑 
块 指 数 (LPIU)、 聚 集 度 指数 (AID) 、 景 观 形状 指数 
(LSI) .香农 多 样 性 指数 (SHDI) 、 香 农 均 匀 度 指数 
(SHEI) 以 及 蔓延 度 指数 (CONTAG)7 个 指数 ( 表 11) 


作为 评价 指标 ,并 利用 Fragstats 软件 计算 指标 值 。 
2.2.2 景观 生态 风险 指数 构建 ”采用 等 间距 的 方式 
将 玻 勒 河流 域 划 分 了 共 1212 个 10 kmx10 km 的 正 
方形 网 格 评价 单元 。 分 别 计算 每 一 个 评价 单元 的 
生态 风险 指数 ,依次 作为 评价 单元 中 心 点 的 景观 生 
态 风 险 , 并 用 克 里 金 法 进行 玖 勒 河流 域 景观 生态 风 
险 的 空间 插值 。 

借鉴 相关 学 者 研究 3, 用 景观 分 离 度 、 景 观 破 
碎 度 以 及 景观 分 维 数 来 建立 景观 干扰 度 指 数 , 分 别 
将 景观 破碎 度 ,分离 度 和 分 维 数 权重 赋值 为 0.5、 
0.3 .0.2 ,进而 可 计算 出 玻 勒 河流 域 景观 生态 风险 指 
数 ,其 公式 及 生态 学 意义 见 表 2。 


表 1 景观 格局 指数 公式 


Tab.1 Formulas of landscape pattern indices 


景观 格局 指数 公式 意义 
斑 块 个 数 (NP) NP=n, nn 为 景观 要 素 i 的 斑 块 指数 ,反映 景观 中 某 一 广 块 类 型 的 斑 块 总 个 数 。 其 值 越 
大 表明 景观 破碎 度 越 高 。 
聚集 度 指数 (ATI) AI=(g,/maxg,)x 100 gi FU maxgi JAE FRA WI DERRE RICZ T ARER TE CS AA 


最 大 斑 块 指数 (LPI) LPI- max(a,a,4,...a,)/A 


景观 形状 指数 (LSI) ise. 
min(£,) 
香农 多 梯 性 指数 (SHDI) — spi yup 


香农 均匀 度 指数 (SHED SHEI --Ë(p xinp) /ln(m) 


接近 0 时 , 斑 块 类 型 破碎 度 最 大 ,凝聚 程度 越 低 ;AI 约 接近 100 时 , 斑 块 类 型 越 
max 为 景观 中 最 大 斑 块 的 面积 (km) ;4 为 景观 总 面积 (km) ;a 为 第 n 个 斑 块 。 
LPI 反 映 最 大 斑 块 在 景观 中 的 优势 比例 。 
瓦 为 第 ;种 土地 利用 类 型 斑 块 边界 的 总 长 度 (km) ,反映 整体 景观 的 形状 复杂 程 
度 。LSI 越 接近 于 1 ,整体 景观 越 简单 ;LSI 越 大 , 则 斑 块 形状 越 不 规则 , 越 离散 。 
P 为 每 一 种 玉 块 类 型 所 占 景 观 总 面积 的 比例 (%);n 为 景观 中 斑 抉 类 型 的 总 
数 。SHDI 反 映 斑 块 类 型 多 样 性 ,多 样 性 指数 越 高 表明 景观 中 斑 块 类 型 越 多 。 

m 为 土地 利用 类 型 数 ;p; 为 每 一 种 斑 块 类 型 所 占 景 观 总 面积 的 比例 (%)。SHEI 


f 值 较 小 时 ,反映 景观 受 一 种 或 儿 种 优势 斑 块 类 型 所 支配 ; 当 SHEI 趋 于 1 时 ,说 
明 景观 中 没有 明显 的 优势 类 型 且 斑 块 类 型 在 景观 中 均匀 分 布 。 
蔓延 度 指数 (CONTAG ) 5» l(p) 7 为 土地 利用 类 型 数 ;n 为 斑 块 数目 ;pj 为 随机 选择 的 2 个 相 邻 栅 格 细胞 是 属于 
contac- 全 全 ” ”i 和 和 j 的 概率 ;CONTAG 反 映 景观 中 不 同类 型 斑 块 的 团聚 程度 或 延展 趋势 。 
21n(m) 


表 2 景观 生态 风险 指数 公式 


Tab.2 Formulas of Landscape ecological risk indices 


指数 名 称 


公式 


景观 生态 风险 指数 (ER,) » T 
景观 生态 风险 指数 (ER mne (Asa) 
k 


izl 


ER SU E des AUS: PE fr k BE SSOUUXUES TECIUC, [EUER A S VA YE UT ROC HE 8A 
仿 值 越 高 ,反之 , 值 越 低 ;4 为 景观 生态 风险 评价 单元 中 第 i 类 景观 的 面积 


(km?) ;4 为 评价 单元 £ 的 总 面积 (km ) ;LL 为 第 i 类 景观 的 生态 损失 指数 ;m 为 


景观 类 型 的 数量 。 
景观 损失 度 指数 (LL) DiC LL 为 景观 损失 度 指数 ;C; 为 破碎 度 指 数 ;U 为 干扰 度 指数 。 
景观 破碎 度 指数 (C)) TES C; 为 破碎 度 指数 ;ni 为 景观 i 的 斑 块 数 ;4 为 景观 i 的 总 面积 (km’)。 
ET 
景观 分 离 度 指 数 (F) jad [n FON OULAT ES BETE TIUS A 为 景观 总 面积 (km) ;4 为 景观 ;的 总 面积 (km2?) ;nn 为 
724, [A 景观 ;的 斑 块 数 。 
景观 分 维 数 (D) 2w(&) DD; 为 景观 分 维 数 ;0; 为 景观 i 的 周 长 (km) ;4 为 景观 ;的 总 面积 (km )。 
4 
D.= US 
景观 干扰 度 指数 (U) U0; 为 景观 干扰 度 指数 ;Ci 为 破 雄 度 指数 ;为 景观 分 离 度 指数 ;D; 为 景观 分 维 


U,=aXC,+bXF,+cXD, 


Jit sa .ce 分 别 为 景观 破碎 度 .分 离 度 和 分 维 数 权重 ,对 其 赋值 为 0.5 .0.3 .0.2。 
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2.2.3 景观 生态 风险 等 级 划分 ”在 ArcGIS 中 ,根据 
自然 断 点 法 将 流域 景观 生态 风险 等 级 划分 低 景观 
生态 风险 较 低 景观 生态 风险 中 景观 生态 风险 , 较 
高 景观 生态 风险 和 高 景观 生态 风险 5 个 等 级 。 
2.2.4 人 为 干扰 度 ”为 分 析 足 勒 河流 域 人 为 干扰 度 
对 景观 生态 风险 的 影响 ,参考 前 人 的 研究 成 果 
结合 研究 区 的 景观 类 型 和 实际 情况 ,对 研究 区 进行 
人 为 干扰 度 计算 ,公式 为 : 


HD--- (1) 


式 中 :HD 为 某 个 网 格 的 人 为 干扰 度 ; UI 为 第 i 类 景 
观 类 型 的 干扰 度 指 数 ;5; 为 第 i 类 景观 类 型 的 面积 
(km?) ;5 为 网 格 单元 的 总 面积 (km ) ;m 为 景观 类 型 
的 数量 。 
2.2.5 地 理 探测 器 ”地理 探 测 咒 具有 因子 探测 e 
互 探测 、 生 态 探测 和 风险 探测 4 种 功能 ,是 由 王 劲 峰 
等 ”提出 的 一 门 新 的 统计 学 方法 , 相 比 其 他 统计 学 
方法 ,地 理 探测 需 具 有 处 理 各 类 别 变量 的 优势 。 目 
前 ,地 理 探测 器 应 用 于 多 学 科 领 域 。 本 文 使 用 前 3 
个 探测 器 。 地 理 探测 顺 的 原理 是 :把 1 个 区 域 划分 
为 若干 个 小 区 域 ,如 果 这 些小 区 域 的 方差 之 和 小 于 
这 个 区 域 整体 的 方差 ,那么 这 个 区 域 具有 差异 性 ; 
如 果 2 个 变量 的 空间 分 布 相同 ,那么 这 2 个 变量 具 
有 一 定 的 关联 性 。 其 中 用 g 值 来 表达 因子 探测 、 空 
间 差 异 和 因子 之 间 的 关系 。 

因子 探测 絮 主 要 探究 了 因 变 量 (7) 的 空间 分 布 


_ hzl 1- SSW (2) 


式 中 :hh 为 自 变 量 的 分 区 序号 ;LL 为 分 区 总 数 ;和 VN, 和 NN 
分 别 为 每 个 分 区 和 区 域内 总 共 的 栅 格 数目 ; on; 和 
o 为 每 个 分 区 方差 和 区 域内 景观 生态 风险 的 方 
差 。SSW 和 SST 分 别 是 层 内 方差 之 和 及 全 区 总 方 
差 。g 的 值 域 为 [0, 1], 值 越 大 说 明 Y 的 空间 差异 越 
明显 。 

交互 作用 探测 主要 在 地 理 探测 器 进行 空间 著 
加 , 旨 在 判定 不 同 因子 之 间 的 交互 力 大 小 ,可 以 冠 
量 表示 交互 后 的 因子 对 景观 生态 风险 的 影响 力 程 
度 。 生 态 探测 则 通过 两 两 因子 进行 比较 ,判断 芷 勒 
河流 域 景观 生态 风险 的 空间 分 布 是 否 具 有 显著 性 
差异 。 


3 结果 与 分 析 


3.1 景观 格局 分 析 

由 表 3 可 知 ,除草 地 景观 以 外 ,其 他 景观 斑 块 个 
数 (NP) 介 于 100~300, 而 草地 的 NP 大 于 1700。 其 
原因 是 草地 分 布 范围 较 广 ,从 平原 到 高 海拔 地 区 均 
出 现 , 具 有 点 - 面 状 相 结合 分 布 现 象 ,破碎 化 程度 
大 ;其 他 景观 要 素面 积 较 小 且 分 布 较 集中 ,未 利用 
地 面积 虽 最 大 ,但 分 布 集中 ,因此 NP 值 较 小 。 最 大 
斑 块 指数 (LPI) 中 ,未 利用 地 LPI 在 70 以 上 ,草地 
LPI 在 2.6 左 右 , 其 余 土 地 利用 LPI 在 0.03~0.44 之 间 


异 质 性 及 其 自 变 量 (X) 在 多 大 程度 上 解释 了 这 种 分 ”波动 ,这 是 因为 未 利用 地 和 草地 为 流域 的 主要 景观 
异 。 其 公式 为 : 类 型 ,面积 较 大 且 分 布 较 广 。 景 观 形状 指数 (LSI) 
表 3 2000—2018 年 疏 勒 河流 域 斑 块 类 型 景观 格局 指数 统计 
Tab.3 Statistics of landscape pattern indices of patch types in the Shule River Basin from 2000 to 2018 

景观 格局 指数 年 份 耕地 林地 草地 水 域 建设 用 地 未 利用 地 
斑 块 个 数 (NP) 2000 127 273 1736 138 85 281 
2010 130 276 1775 146 84 293 
2018 126 261 1770 159 164 331 

最 大 斑 块 指数 (LPI) 2000 0.3014 0.0916 2.6584 0.1058 0.0373 78.1309 

2010 0.3370 0.0916 2.6602 0.1058 0.0373 77.7646 

2018 0.4360 0.0810 2.7369 0.0997 0.0391 77.3710 

景观 形状 指数 (LSI) 2000 16.9571 17.549 58.4394 14.5510 10.4194 30.4054 

2010 17.6835 18.0000 58.9516 14.7000 10.3548 30.7383 

2018 17.4118 17.6863 58.4410 15.2500 12.9730 31.1059 

聚集 度 指数 (AT) 2000 52.1012 30.1902 59.7878 41.2909 31.6159 90.1088 

2010 56.1398 29.1667 59.5156 40.9483 32.4009 89.9689 

2018 59.8792 30.5306 59.6768 41.0032 31.5301 89.8246 
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中 ,草地 > 未 利用 地 > 林地 > 耕地 > 水 域 > 建设 用 地 ,由 
于 草地 生长 范围 广 ,多 分 布 在 未 利用 地 周围 ,而 其 
他 景观 要 素 ( 如 建设 用 地 和 水 域 ) 存 在 人 为 干涉 进 
行规 划 等 限制 手段 ,因此 草地 、 未 利用 地 斑 块 形状 
不 规则 ,建设 用 地 等 斑 块 形状 较 规 则 。 从 聚集 度 指 
数 (AIT) 看 出 ,未 利用 地 AI 在 90 左 右 , 表 明 未 利用 地 
斑 块 聚集 程度 高 ,其 余 景观 要 素 AI 介 于 29.17~ 
59.88。 

由 表 4 可 知 , 玻 勒 河流 域 景观 类 型 蔓延 度 指数 
在 2000 一 2018 年 时 下 降 趋势 ,说 明 这 些 年 内 , 随 着 
人 口 不 断 增 加 以 及 土地 开 暴 . 扎 荒 等 土地 利用 变 
化 ,景观 破碎 程度 不 断 增加 ;而 香农 多 样 性 指数 和 
香农 均匀 度 指 数 在 2000 一 2018 年 呈 缓 慢 上 升 趋势 ， 
但 值 偏 小 ,说 明 斑 块 分 布 仍 不 均匀 。 


表 4 2000—2018 FRANTA ARME 
景观 格局 指数 统计 
Tab.4 Statistics of landscape pattern indices of landscape 
types in the Shule River Basin from 2000 to 2018 


PUN 蔓延 度 指数 香农 多 样 性 香农 均匀 度 
Di (CONTAG) ”指数 (SHDI1) EP (SHED) 
2000 69.9751 0.6147 0.3431 
2010 69.3664 0.6284 0.3507 
2018 68.6688 0.6434 0.3591 


3.2. 景观 生态 风险 分 析 

从 空间 上 看 ,2000 一 2018 年 琉 勒 河流 域 景观 生 
态 风 险 等 级 (图 2) ,分布 规律 趋势 一 致 ,高 景观 生态 
风险 等 级 主要 分 布 在 流域 北部 , 低 景 观 生态 风险 等 
级 分 布 在 流域 南部 。 原 因 是 流域 北部 景观 类 型 主 
要 为 未 利用 地 ,景观 类 型 单一 旦 对 外 界 的 抵抗 能 
差 , 从 而 生态 环境 相对 脆弱 ,而 流域 南部 由 于 靠近 
祁连山 脉 ,水 资源 丰富 、 景 观 类 型 多 样 ,主要 有 草 
地 水域 及 林地 ,草地 及 林地 水 源 涵 养 能 力 强 ,南部 
海拔 较 高 ,人 类 干扰 少 , 因 此 景观 生态 风险 等 级 低 。 

由 表 5 可 知 ,2000 一 2010 年 整体 景观 生态 风险 
等 级 变化 不 大 ,主要 在 较 高 和 高 景观 生态 风险 等 级 
之 间 转 换 。 这 一 时 段 内 ,未 利用 地 和 林地 面积 减 
少 ,耕地 草地、 建设 用 地 面积 增加 ,其 中 ,耕地 面积 
增加 351.51 km ,这 是 由 于 移民 工程 导致 的 耕地 扩 
张 , 人 类 对 土地 开垦 现象 严重 ,导致 景观 生态 风险 
向 高 值 转移 。 

2010 一 2018 年 , 较 高 和 高 景观 生态 风险 等 级 面 
积 减 少 ,其 中 较 高 景观 生态 风险 等 级 面积 减少 
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图 2 2000 一 2018 年 牙 勒 河流 域 景 观 生态 风险 等 级 分 布 
Fig. 2 Distributions oflandscape ecological risk levels in the 
Shule River Basin from 2000 to 2018 
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表 5 2000 一 2018 年 琉 勒 河流 域 景观 生态 风险 等 级 面积 及 比例 统计 
Tab.5 Statistics of the area and proportion of landscape ecological risk grades in the Shule River Basin from 2000 to 2018 


景观 生态 2000 一 2010 年 2010 一 2018 年 2000 一 2018 年 

风险 等 级 面积 /km 比例 /% 面积 km 比例 /% 面积 /km 比例 /% 
低 -23.51 -0.02 1154.71 1.03 1131.19 1.01 
较 低 -37.16 -0.03 1071.26 0.95 1034.10 0.92 
m -5.03 -0.01 2231.55 1.99 2226.52 1.98 
较 高 -482.94 -0.43 -3230.81 -2.88 -3713.75 -3.30 
高 545.51 0.49 -1228.21 -1.09 -682.71 -0.61 
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2.8896 , 主要 在 玉 门 市 西南 部 ,高 景观 生态 风险 等 级 
面积 减少 1.09% ,主要 在 肃 北 ( 南 ) 中 部 一 带 , 而 低 、 
较 低 及 中 景观 生态 风险 等 级 比例 在 瓜 州 县 中 部 及 
肃 北 ( 南 ) 中 部 等 地 开始 上 升 。 由 于 流域 南部 主要 
为 高 海拔 地 区 ,水 域 . 草 地 、 林 地 景观 面积 增加 ,说 
明 该 区 域 景 观 生态 风险 有 所 降低 ,其 主要 原因 是 加 
强 了 对 祁连山 生态 环境 保护 ,下 游 地 区 开始 实施 节 
水 战略 ,如 根据 《敦煌 水 资源 合理 利用 与 生态 保护 
综合 规划 》, 敦 煌 实施 了 灌区 节 水 改造 种植 结构 调 
整 等 措施 ,遏制 了 敦煌 盆地 地 下 水 下 降 趋势 ,使 得 
高 景观 生态 风险 等 级 逐渐 向 中 、 较 低 甚至 低 等 级 转 
移 。 瓜 州 县 中 部 为 蕊 勒 河 谷地 ,景观 生态 风险 等 级 
降低 。 近 18 a, 较 高 景观 生态 风险 等 级 面积 下 降 
3.30% ,高 景观 生态 风险 面积 下 降 0.61% , 低 、 较 低 和 
中 景观 生态 风险 等 级 面积 分 别 上 升 1.01% .0.92% 和 
1.98%, 

2000—2010 4E. , E 2i Je b T A RESA 
险 等 级 区 的 面积 有 所 增加 ,增加 了 545.51 km’, 占 流 
域 比 为 0.49%, 低 、 较 低 、 中 等 和 较 高 景观 生态 风险 
等 级 区 的 面积 均 有 所 下 降 ,分 别 减少 23.51 km’, 
37.16 km 、5.03 km 和 482.94 km’, 占 流域 面积 的 比 
例 分 别 减少 了 0.02% 、0.03%、0.01% 和 0.43%。 从 空 
间 分 布 (图 3) 来 看 ,景观 生态 风险 增加 的 区 域 零星 
分 布 在 流域 各 处 ,景观 生态 风险 等 级 下 降 的 区 域 则 
主要 集中 在 敦煌 市 和 玉 门 市 ,由 《中 国 县 域 统计 年 
鉴 ) 统 计 得 出 ,敦煌 市 人 口 由 2000 年 的 1.7x10; 人 减 
少 到 2010 年 的 1.4x10 人, 玉 门市 由 1.9x10 人 减少 
为 1.5x10; 人 ,由 于 人 口 减少 导致 的 人 类 干扰 下 降 ， 
加 上 政府 对 玉 门 资源 开采 的 控制 ,导致 景观 生态 风 
险 等 级 下 降 。 

2010 一 2018 年 ,流域 低 景观 生态 风险 面积 增加 
1154.71 km' , 较 低 景观 生态 风险 面积 增加 1071.26 
km ,中 景观 生态 风险 等 级 面积 增加 2231.55 km , 较 


高 和 高 景观 生态 风险 等 级 面积 共 减 少 4459.02 km’, 
整个 流域 景观 生态 风险 等 级 逐渐 下 降 , 且 降低 区 域 
相对 集中 ,具体 分 布 在 肃 北 ( 北 ) 中 部 及 西部 .敦煌 
市 中 部 、 瓜 州 县 东部 及 玉 门 市 西部 。 在 此 期 间 , 玉 
门市 规模 以 上 工业 企业 数 由 50 个 减少 为 4 个 ,部 
煌 市 由 46 个 减少 为 23 个 , 肃 北 由 24 个 减少 为 14 
个 , 瓜 州 县 由 32 个 增加 到 41 个 ,但 除 瓜 州 县 人 口 减 
少 1x10* 人 外 ,其 他 县 市 人 口 均 稳 定 , 人 类 干扰 强度 
减弱 ,流域 整体 景观 生态 风险 等 级 下 降 。 

总 体 来 说 ,2000 一 2018 年 流域 低 景观 生态 风险 
面积 增加 1131.19 km? , 较 低 景观 生态 风险 面积 增加 
1034.10 km’, 中 景观 生态 风险 等 级 面积 增加 
2226.52. km , 较 高 和 高 景观 生态 风险 等 级 面积 共 减 
少 4396.46 km’。 

3.3 景观 生态 风险 空间 差异 驱动 因子 分 析 

因子 探测 器 既 能 探测 玻 勒 河流 域 景 观 生态 风 
险 的 空间 分 异性 ,也 能 探测 影响 因子 在 多 大 程度 上 
解释 景观 生态 风险 空间 分 异性 。 为 方便 计算 ,Xi、 
Xi Xi Xi XXX 和 Xs 分 别 代 表 高 程 .坡度 降水 、 
NDVI, 最 高 气温 、 最 低 气 温 .人 口 密度 和 人 为 干 
扰 度 。 

由 图 4 可 知 ,2000 年 决定 力 g 值 从 大 到 小 依次 
为 Xs>X4>Xs>Xs>Xi>Xe>X>X ,人 为 干扰 度 在 影响 景 
观 生态 风险 变化 的 贡献 量 中 排 在 首位 ,说 明 人 为 干 
扰 度 是 琉 勒 河流 域 景 观 生 态 风险 变化 的 主要 控制 
因素 ,第 二 主导 因素 是 NDVL, 贡献 率 为 0.2566 ,其 次 
主导 因素 为 高 程 . 降 水 .最 高 和 最 低 气 温 , 人 口 密度 
是 解释 力 最 弱 的 一 个 因素 ,贡献 率 仅 为 0.0142。 
2010 年 的 4 值 排列 变化 与 2000 年 一 致 。2018 年 与 
前 2 个 时 期 的 排列 顺序 稍 有 不 同 ,X, 由 前 2 个 时 期 
的 第 五 排 到 第 三 ,X; 则 由 第 三 降 到 第 五 。3 个 时 期 g 
值 对 比 可 知 ,不 同 因子 在 不 同年 份 变化 不 同 ,X~Xs 
的 g 值 在 2000 一 2018 年 均 呈 现 “V” 型 变化 , 即 6 类 因 
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(a) 2000—20104E 


景观 生态 风险 空间 变化 


国 景观 生态 风险 等 级 降低 0 70 km 
O 景观 生态 风险 等 级 不 变 — 
ES. 景观 生态 风险 等 级 升 高 

(b) 2010 一 2018 年 


Xz 


景观 生态 风险 空间 变化 

国 景 观 生 态 风险 等 级 降低 0 70 km 
CE 景观 生态 风险 等 级 不 变 [| 
国 四 景观 生态 风险 等 级 升 高 


(c) 2000 一 2018 年 


景观 生态 风险 空间 变化 

国 景观 生态 风险 等 级 降低 
上 喇 景观 生态 风险 等 级 不 变 
ENS 景观 生态 风险 等 级 升 高 


图 3 2000 一 2018 年 芷 勒 河流 域 景观 生态 风险 空间 变化 
Fig.3 Spatial changes of landscape ecological risk in the 
Shule River Basin from 2000 to 2018 
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DE X 6 A X Xo、 各 和 XX 分 别 代 表 高 程 坡 度 、 降 水 、NDVI、 最 
高 气温 、 最 低 气 温 、 人 口 密度 和 人 为 干扰 度 。 下 同 。 
图 4 2000 一 2018 年 巩 勒 河流 域 景观 生态 风险 因子 解释 力 
Fig.4 Explanatory power of landscape ecological risk 
factors in the Shule River Basin from 2000 to 2018 


子 在 2010 年 g 值 最 小 ,2018 年 g 值 最 高 。X~X 则 呈 
现 “ 和 ”型 变化 ,2010 年 g 值 最 高 。 整 体 上 ,2000 一 
2010 年 各 因子 决定 力 9 值 接近 ,在 2010 一 2018 年 4 
值 变 化 波动 大 ,这 与 3 个 时 期 景观 生态 风险 等 级 分 
布 图 变化 一 致 。 

玻 勒 河流 域 景观 生态 风险 空间 分 异 的 驱动 因 
子 交 互 探测 结果 显示 ( 表 6) ,所 选 因子 中 任意 2 个 因 
子 的 交互 作用 均 大 于 单个 因子 ,交互 类 型 主要 为 双 
因子 增强 和 非 线性 增强 2 种 ,说明 玻 勒 河流 域 景观 
生态 风险 空间 分 异 的 结果 不 是 由 单一 影响 因子 造 
成 的 ,而 是 由 多 影响 因素 共同 作用 的 结果 。 其 中 ， 
2000 年 夭 与 和 交互 作用 的 影响 最 强 ,g 值 最 高 为 
0.6848 ; 交互 影响 程度 达到 0.6 以 上 的 还 有 :Xs 站 Xi、 
XX 、X; 和 Xe。 其 他 因子 交互 9 值 虽然 小 于 0.6, 但 
也 显示 了 双 因 子 比 单 因子 对 于 景观 生态 风险 影响 
程度 更 高 。2010 年 ,X 与 和 交互 作用 最 强 ,g 值 为 
0.7178, 其 他 因子 与 夭 交 互 时 ,g 值 均 在 0.6 以 上 。 
2018 年 ,XX 与 和 交互 时 ,g 值 最 高 , 值 为 0.5757 ,自然 
因子 之 间 交 互 作用 增强 ,但 贡献 率 均 在 0.45 以 下 ， 
只 有 自然 因子 和 人 类 活动 因子 共同 作用 时 ,对 景观 
生态 风险 的 解释 更 强 ,XX 作为 主导 因子 影响 琉 勒 河 
流域 景观 生态 风险 ,同时 与 气候 及 地 表 因 素 的 共同 
作用 下 ,对 景观 生态 风险 分 布 影响 最 大 。 

玖 勒 河流 域 景观 生态 风险 空间 分 析 驱 动因 子 
之 间 的 交互 作用 均 大 于 单 因 子 对 景观 生态 风险 的 
作用 程度 ,不同 因子 之 间 的 作用 所 形成 的 增强 效应 
共同 影响 了 流域 景观 生态 风险 在 空间 上 分 异 效 
果 。 其 中 ,3 个 时 期 人 为 干扰 度 与 其 他 因子 的 交互 
作用 均 最 高 ,明显 表现 出 多 因子 作用 增强 的 效果 。 
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孙 丽 车 等 :基于 地 理 探 测 咒 模型 的 下 勒 河流 域 景观 生态 


风险 评价 及 驱动 因素 分 析 


表 6 2000 一 2018 年 琉 勒 河流 域 景观 生态 风险 驱动 因子 交互 影响 


Tab.6 Interaction of driving factors of landscape ecological risk in the Shule River Basin from 2000 to 2018 


2000 年 2010 年 2018 年 

交互 项 KHE ”交互 项 AH 交互 项 KEE ”交互 项 交互 值 交互 项 ZHE ”交互 项 ”交互 值 
XNX, 0.2474 XMXs 0.3824 XNX, 0.2421 XMXs — 0.3576 XNX, | 02554 XNX;, — 02962 
XXX, 0.3492 X4X, — 0.324T XNX; — 0.3327 X4X, — 0322V XX, 0.3163 XNXs — 0.2993 
XX, . 0.3760 XMX — 02074 XX, 0.3692 X4X, — 02072 XX, | 04459 XX, | 0262l 
XX, | 02479 XMNXs —— 0.6848 XX. 02334 XMXs — 0.7024 XX. 0.2909 XX, — 0.5011 
XX, | 02364 XMXs | 0.3960 XX, 02327 XNMXs | 0.3769 XX, 0.3028 XX. — 04315 
XNX, 020902 . XX, 079r XX, | 0.2082 XX, | 0.3661 XX, 02713 XANAX. 04248 
XNMXs ^ 0.6456 XNMX | 02770 XNXs 0.6934 X(XX, | 02974 XImXs 0.5239 XX, | 02735 
XX. — 02423 XX. —— 0.6449 XX. — 02464 XimX。 — 0.6491 XX. — 02613 XX. — 0.5757 
XX, — 0.3166 X4X, | 0.2684 XX, —— 0.3249 X4X, — 02513 XX, — 0.276 XNXs — 02825 
XX; | 02680 XMX — 02319 XX. | 02535 XMX — 02126 XNX; — 0.2499 XNX, — 0.2686 
XMXs | 02352 XX:s —— 0.6845 XX, 02227 XANAX; — 07178 XX. — 02489 XNX: — 0.5178 
XX, — 0.136 — X41X; — 0.1504 XX, — 0.1064 XANAX, — 0.1250 XX, 01663 XAMX; — 02014 
XMX: —— 0.6507 XANAX: 0.6665 XX, — 0.6794 XANAX: —— 0.6972 XX. —— 0.4608 XANXs — 0.5333 
XMX, — 04148 XMX: —— 0.5998 XX, — 04221 X4X, — 0.6337 XX, — 04325 XAMXs — 04483 

B: Xn Xa Xa Xa XXX 和 Xs 分 别 代表 高 程 . 坡 度 、 降 水 NDVI` 最 高 气温 、 最 低 气 温 `、 人 口 密度 和 人 为 干扰 度 ;* 代 表 双 因子 增强 。 

从 交互 结果 来 看 ,人 类 活动 因子 与 自然 因子 作用 明 — 均 存在 显著 性 差异 ,不 存在 显著 性 差异 主要 表现 在 


显 强 于 自然 因子 之 间 交 互 。 
生态 探测 反映 每 2 个 因子 对 景观 生态 风险 空 


分 布 的 影响 是 否 具 有 显著 差异 ,可 以 进 一 fee 
导 的 影响 因子 ,并 评价 其 作用 机 理 的 差异 性 。2000 


年 生态 探测 图 显示 XmX X; Xs X; X XX XN 
X5 X, 1X5 为 无 显著 差异 ;2010 年 生态 探测 图 显示 


自然 因子 内 部 。 
4 讨论 


4.1 景观 生态 风险 时 空 分 异 特征 
ifc) 河流 域 景 观 生 态 风 险 空 间 分 布 北 高 南 低 ， 


Xi iso IX eoe Xa X OG SIX S 
差异 ;2018 年 生态 探测 图 显示 Xi 站 Xs XX Xi Xs 
KNX 、Xs 门 Xs、Xs 门 Xs 、X 门 Xs; 为 无 显著 差异 (图 5)。 
从 3 个 时 期 的 显著 性 差异 可 以 看 出 ,X; 与 其 他 因子 


"m 


(a) 2000 年 


(b) 20104] 


1 
X, X, X, X, X; X, X; X, X, X, X 
驱动 因子 


与 潘 竟 虎 等 汪 人 研究 结果 一 致 ,流域 从 北向 南 ,海拔 
逐渐 升 高 ,温度 逐渐 降低 ,降水 逐渐 增加 ,景观 类 型 
从 未 利用 地 .荒漠 耕地、 林地 过 渡 到 草地 冰川 , 景 
观 生态 风险 随 之 降低 。2000 一 2010 年 ,整体 景观 生 


(c) 2018 年 


X, X; X; X; Xs X; XQ X. XQ X. X2 X; XQ 
驱动 因子 


驱动 因子 


注 :Y 为 存在 显著 性 差异 ;N 为 不 存在 显著 性 差异 。 
图 5 2000 一 2018 年 玻 勒 河流 域 景观 生态 


风险 驱动 因子 显著 性 差异 


Fig. 5 Significant differences in driving factors of landscape ecological risk in the Shule River Basin from 2000 to 2018 
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态 风险 等 级 变化 不 大 ,主要 在 较 高 风险 等 级 和 高 景 
观 生 态 风 险 等 级 之 间 转 换 , 景 观 生态 风险 增加 的 区 
域 零 星 散 状 分 布 ,景观 生态 风险 等 级 下 降 的 区 域 则 
主要 集中 在 敦煌 市 和 玉 门 市 。2010 一 2018 年 , 较 高 
和 高 景观 生态 风险 等 级 面积 减少 3.97% ,表明 流域 
生态 环境 在 整体 上 呈现 好 转 趋势 。 康 紫 蔽 等 “在 
同样 处 于 干旱 内 陆 区 的 玛 纳 斯 河流 域 研究 发 现 类 
似 结果 。 造 成 景观 生态 风险 高 或 者 较 高 主要 是 受 
未 利用 地 和 城建 用 地 影响 ,而 景观 生态 风险 等 级 
较 低 和 低 的 景观 类 型 主要 受 林地 .草地 及 园地 影 
响 ,这 与 在 博 斯 腾 湖 流域 所 、 白 龙 江 流域 ”红河 
流域 全 和 梨 湖 流域 “等 不 同 生态 条 件 的 研究 结果 
一 致 ;而 与 细 河 流域 ”不同 ,其 主要 是 该 流域 北部 
工矿 建设 用 地 面积 相对 集中 ,受到 人 类 干扰 后 的 损 
失 度 低 , 因 此 景观 生态 风险 低 。 划 分 景观 生态 风险 
网 格 ,构建 景观 生态 风险 指数 揭示 了 玻 勒 河流 域 景 
观 生态 风险 时 空 特征 ,与 研究 区 生态 环境 相 吻 合 ， 
说 明 在 较 短 时 间 景观 生态 风险 能 够 反映 人 类 活动 
对 其 影响 ,但 流域 是 一 个 复杂 的 综合 体 , 虽 通 过 景 
观 生态 风险 网 格 客观 反映 景观 生态 风险 ,但 网 格 划 
分 是 以 等 间距 划分 ,对 考虑 流域 内 干 湿 、 农 牧 交 葵 
区 “等 变化 缺乏 考察 , 须 在 日 后 研究 中 进行 完善 。 
4.2 驱动 因素 

驱动 因子 分 析 表 明 人 为 干扰 度 是 影响 景观 生 
态 风 险 空间 分 异 的 主要 控制 因子 ,人 口 密度 影响 最 
弱 ,与 买 世 军 光 对 喀斯特 山区 流域 综合 生态 风险 红 
动力 分 析 结 果 相 同 。 因 子 交 互 作用 结果 均 为 双 因 
子 增强 和 非 线 性 增强 ,说 明 因子 间 两 两 交互 作用 比 
单 因子 作用 对 景观 生态 风险 的 影响 更 强 ,尤其 表现 
在 人 为 干扰 度 与 其 他 因子 交互 。 人 为 干扰 度 与 其 
他 因子 均 存 在 显著 性 差异 ,不 存在 显著 性 差异 主要 
表现 在 自然 因子 。 郑 续 等 中 研究 该 区 域 产 水 量 时 
空 特征 时 也 发 现 人 为 干扰 强度 与 其 他 因子 之 间 的 
交互 作用 显著 相关 。 由 此 可 见 , 对 于 下 勒 河流 域 景 
观 生态 风险 调控 来 说 主要 是 控制 人 类 干扰 度 。 鉴 
于 人 为 活动 因子 包括 一 些 政策 性 因素 等 由 于 数据 
获取 和 定量 化 具有 一 定 难 度 ,因此 本 研究 并 未 过 多 
考虑 ,实际 上 不 同 区 域 经 济 调控 、 生 态 保护 政策 的 
差异 对 景观 生态 风险 的 影响 是 巨大 的 ,这 也 是 未 来 
研究 的 重点 问题 之 一 。 


5 结论 


< 


本 文 以 玻 勒 河流 域 为 研究 对 象 ,构建 景观 生态 
风险 指数 和 利用 地 理 探 测 器 分 析 影 响 玻 勒 河流 域 
景观 生态 风险 的 驱动 因素 ,得 出 以 下 主要 结论 : 

(1) 通过 对 2000 一 2018 年 芍 勒 河流 域 王 块 类 
型 景观 格局 指数 的 分 析 ,可 以 得 到 草地 的 斑 块 个 数 
(CNP) 和 景观 形状 指数 (LST) 最 高 ,未 利用 地 的 最 大 
斑 块 指数 (LPI) 和 聚集 度 指数 (AI) 最 高 ,主要 原因 是 
草地 和 未 利用 地 作为 流域 的 主要 景观 类 型 ,面积 较 
大 且 分 布 较 广 。2000 一 2018 Fi 1$] Ai Je c XUL E 
延 度 指数 (CONTAG) 呈 下 降 趋势 ,而 多 样 性 指数 
(SHDIT) 和 均匀 度 指 数 (SHET) 呈 缓慢 上 升 趋势 ,说 明 
景观 破碎 程度 不 断 上 升 。 

(2) 2000 一 2018 年 芍 勒 河流 域 景观 生态 风险 空 
间 呈 现 北 高 南 低 分 布 特征 。 其 中 , 较 高 景观 生态 风 
念 等 级 下 降 3.30%, 高 景观 生态 风险 下 降 0.61%, 低 、 
较 低 和 中 景观 生态 风险 等 级 分 别 上 升 1.01%、0.92% 
和 1.98%。 

(3) 人 为 干扰 度 是 影响 景观 生态 风险 空间 分 布 
的 主要 因子 ,NDVI 次 之 ,而 人 口 密度 影响 最 小 。 影 
响 因 子 对 景观 生态 风险 的 影响 不 是 独立 的 ,而 是 因 
子 共 同 作用 的 结果 。 不 同 影响 因子 之 间 存 在 交互 
作用 ,是 任意 2 种 影响 因子 的 作用 值 都 大 于 任意 一 
类 单 因子 的 值 ,其 中 人 为 干扰 因子 与 其 他 因子 交互 
作用 显著 ,对 景观 生态 风险 的 影响 均 为 双 因 子 增强 
或 非 线 性 增强 。 
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Abstract: The Shule River Basin is a typical arid inland river basin with an extremely fragile ecological 
environment. It is also one of the most important ecological barriers of northwest China. This study analyzes the 
spatial and temporal variation characteristics of the landscape ecological risk in the Shule River Basin with land 
use data from 2000 to 2018 using Fragstats software. First, the results show that the main landscape types in the 
Shule River Basin are unused land and grassland with high number of patches (NP), landscape shape index (LSI), 
largest patch index (LPI), and aggregation index (AI) values. The contagion (CONTAG) value dramatically 
declined, while the Shannon’ s diversity index (SHDI) and Shannon’ s evenness index (SHEI) values slowly 
increased during the period from 2000 to 2018. This indicates that the landscape fragmentation is serious in the 
basin. Second, the landscape ecological risk in the Shule River Basin in the north was higher than that in the 
south, and a gradual downward trend was observed from 2000 to 2018, especially for the high-risk area. Third, 
human disturbance is the main factor affecting the spatial distribution of the landscape ecological risk, followed 
by NDVI and population density. The impact of the landscape ecological risk was double factor or nonlinear 
enhancement. No significant differences were found between natural and natural factors, while significant 
differences were determined between natural and anthropogenic factors. 
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